Extragalactic Absorption of High Energy Gamma-Rays 
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The pair-production absorption of high-energy 7-rays by intergalactic low-energy photons is 
expected to produce a high-energy cutoff in the spectra of sources which is a sensitive function of 
redshift. We first discuss the expected absorption coefficient as a function of energy and redshift 
derived by Stecker and De Jager by making use of a new empirically based calculation of the spectral 
energy distribution of the intergalactic infrared radiation field as given by Malkan and Stecker. We 
then discuss the fact that new data on the high energy 7-ray source Mrk 501 appear to show the 
amount of intergalactic absorption predicted. The implications of this new HEGRA data, should 
they be confirmed, are significant for the astrophysics of this source, implying that (1) there is 
no significant intrinsic absorption inside the source, and (2) the physics of the emission spectrum 
produces a power-law to energies above 20 TeV. As a further test for intergalactic absorption, we 
give a predicted spectrum, with absorption included, for PKS 2155-304. This XBL lies at a redshift 
of 0.12, the highest redshift source yet observed at an energy above 0.3 TeV. We also discuss the 
determination of the 7-ray opacity at higher redshifts (out to z = 3), following the treatment of 
Salamon and Stecker. 



I. INTRODUCTION 



Very high energy 7-ray beams from blazars can be used to measure the intergalactic infrared radiation field, since 
pair-production interactions of 7-rays with intergalactic IR photons will attenuate the high-energy ends of blazar 
spectra [jlj. In recent years, this concept has been used successfully to place upper limits on the the intergalactic IR 
field (IIRF) Q - Determining the (IIRF), in turn, allows us to model the evolution of the galaxies which produce 
it. As energy thresholds are lowered in both existing and planned ground-based air Cherenkov light detectors 0, 
cutoffs in the 7-ray spectra of more distant blazars are expected, owing to extinction by the IIRF. These can be used 
Oh; to explore the redshift dependence of the IIRF M , M . 

There are now 66 "grazars" (7-ray blazars) which have been detected by the EGRET team JlOfl . These sources, 
optically violent variable quasars and BL Lac objects, have been detected out to a redshift greater that 2. Of all of 
the blazars detected by EGRET, only the low-redshift BL Lac, Mrk 421 (z = 0.031), has been seen by the Whipple 
telescope |jTjj| . The fact that the Whipple team did not detect the much brighter EGRET source, 3C279, at TeV 
energies |fl2"| , fl3]| is consistent with the predictions of a cutoff for a source at its much higher redshift of 0.54 [Q. So 
too are the further detections of three other close BL Lacs (z < 0.12), viz., Mrk 501 (z = 0.034) @, 1ES2344+514 
(z = 0.044) ||, and PKS 2155-304 (z = 0.117) @ which were too faint at GeV energies to be seen by EGREl^. 

The formulae relevant to absorption calculations involving pair-production are given and discussed in Ref. [[l]. For 
7-rays in the TeV energy range, the pair-production cross section is maximized when the soft photon energy is in the 
infrared range: 

A(£ 7 ) ~ A e — = 2AE 1 . T eV fJ-m (1) 

where \ e — hj (m e c) is the Compton wavelength of the electron. For a 1 TeV 7-ray, this corresponds to a soft photon 
having a wavelength near the K-band (2.2y«m). (Pair-production interactions actually take place with photons over a 
range of wavelengths around the optimal value as determined by the energy dependence of the cross section; see eq. 
(3).) If the emission spectrum of an extragalactic source extends beyond 20 TeV, then the extragalactic infrared field 



X PKS 2155-304 was seen in one observing period by EGRET as reported in the Third EGRET Catalogue JlC 
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should cut off the observed spectrum between ~ 20 GeV and ~ 20 TeV, depending on the redshift of the source 0, 



II. ABSORPTION OF GAMMA-RAYS AT LOW REDSHIFTS 

Stecker and De Jager |l7]] (hereafter SD98) have recalculated the absorption coefficient of intergalactic space using 
a new, empirically based calculation of the spectral energy distribution (SED) of intergalactic low energy photons by 
Malkan and Stecker |Q (hereafter MS98) obtained by integrating luminosity dependent infrared spectra of galaxies 
over their luminosity and redshift distributions. After giving their results on the 7-ray optical depth as a function 
of energy and redshift out to a redshift of 0.3, SD98 applied their calculations by comparing their results with the 
spectral data on Mrk 421 @ and spectral data on Mrk 501 @. 

SD98 make the reasonable simplifying assumption that the IIRF is basically in place at a redshifts < 0.3, having been 
produced primarily at higher redshifts |21jj . Therefore SD98 limited their calculations to z < 0.3. (The calculation of 
7-ray opacity at higher redshifts ||, || will be discussed in the next section.) 

SD98 assumed for the IIRF, two of the SEDs given in MS98 Jl8|. Their upper curve now appears to be in better 
agreement with lower limits from galaxy counts, with Keck telescope, HST. NICMOS, ISO and SCUBA studies of 
galaxies at high redshifts (Ref . [plj and references therein) and with COBE data p2| - |27| (see Figure P . 




FIG. 1. The upper infrared SED from Malkan and Stecker compared with observational data and other constraints (courtesy 
O.C. De Jager). 

The results of MS98 are also in agreement with upper limits obtained from TeV 7-ray studies j|| - ||. This 
agreement is illustrated in Figure [l] which shows the upper SED curve from MS98 in comparison with various data 
and limits. 

The SD98 results for the absorption coefficient as a function of energy do not differ dramatically from those obtained 
previously ^8|, p9[| ; however, they are more reliable because they are based on the empirically derived IIRF given 
by MS98, whereas all previous calculations of TeV 7-ray absorption were based on theoretical modeling of the IIRF. 
The MS98 calculation was based on data from nearly 3000 IRAS galaxies. These data included (1) the luminosity 
dependent infrared SEDs of galaxies, (2) the 60/im luminosity function of galaxies and, (3) the redshift distribution 
of galaxies. 
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FIG. 2. Optical depth versus energy for 7-rays originating at various redshifts obtained using the SEDs corresponding to the 
lower IIRF (solid lines) and higher IIRF (dashed lines) levels shown in a Figure taken from SD98. As discussed in the text, the 
higher IIRF curves (dashed lines) are more in line with recent data. 

The advantage of using empirical data to construct the SED of the IIRF, as done in MS98, is particularly indicated 
in the mid-ir range where galaxy observations indicate more flux from warm dust in galaxies than that taken account 
of in more theoretically oriented models. As a consequence, the mid-IR "valley" between the cold dust peak in the 
far IR and cool star peak in the near IR is partially filled in (see Figure [[]). For a source at low redshift, it follows 
from cq. (1) that 7-rays of energy ~ 20 TeV will be absorbed preferentially by photons in the wavelength range of 
this "valley", i.e., near 50 /im. In this range, significant lower limits now exist which are near the predicted IIRF flux 
(see Figure||). 

In fact, the observed flaring spectrum of Mrk 501 has been newly extended to an energy of 24 TeV by observations 
of the HEGRA group [[30|. The new HEGRA data are well fitted by an E~ 2 source spectrum steepened at energies 
above a few TeV by intergalactic absorption with the optical depth calculated by SD98 |5l| . Figure ^, taken from Ref. 
pl| , clearly shows this. The philosophy behind Ref. is that the existing lower limits on the mid-ir background 
flux predict a minimum expected absorption. The derived unabsorbed source spectrum then tells us (1) that there is 
negligible intrinsic absorption in the source, and (2) the physics of the emission mechanism should give a power-law 
spectrum with a spectral index of ~2 up to an energy of at least ~ 20 TeV. 

Consider the source PKS 2155-304, an XBL located at a moderate redshift of 0.117, which has been reported by 
the Durham group to have a flux above 0.3 TeV of ~ 4 x 10~ n cm~ 2 s _1 |16| . We predict that this source should 
have its spectrum steepened by ~ 1 in its spectral index between ^0.3 and ~ 3 TeV and should show an absorption 
turnover above ~ 6 TeV as shown in Figure || Observations of the spectrum of this source should provide a further 
test for intergalactic absorption. 

III. ABSORPTION OF GAMMA-RAYS AT HIGH REDSHIFTS 

In order to calculate high-redshift absorption properly, it is necessary to determine the spectral distribution of 
the intergalactic low energy photon background radiation as a function of redshift as realistically as possible out 
to frequncies beyond the Lyman limit. This calculation, in turn, requires observationally based information on the 
evolution of the spectral energy distributions (SEDs) of IR through UV starlight from galaxies, particularly at high 
redshifts. 

Salamon and Stecker || (hereafter SS98) have calculated the 7-ray opacity as a function of both energy and redshift 
for redshifts as high as 3 by taking account of the evolution of both the stellar population spectra and emissivity 
of galaxies with redshift. In order to accomplish this, they adopted the recent analysis of Fall, et al. [53 and also 
included the effects of metallicity evolution on galactic stellar population spectra. They then gave predicted 7-ray 
spectra for selected blazars and extend our calculations of the extragalactic 7-ray background from blazars to an 
energy of 500 GeV with absorption effects included. 
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FIG. 3. The bottom curve and points show the new HEGRA data on Mrk 501 in the flaring state fl3(i|| ; the upper line and 
points show the intrinsic spectrum ol the source with the effect of the predicted extragalactic absorption removed |3jj|. 
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FIG. 4. Predicted differential absorbed spectrum, for PKS 2155-304 (solid line) assuming an E 2 differential source spectrum 
(dashed line) normalized to the integral flux given in Ref. JlfJ (see text). 

Fall, et al. |32| have devised a method for calculating stellar emissivity which bypasses the uncertainties associated 
with estimates of poorly defined luminosity distributions of evolving galaxies. The core idea of their approach is to 
relate the star formation rate directly to the evolution of the neutral gas density in damped Ly a systems, and then 
to use stellar population synthesis models to estimate the mean co- moving stellar emissivity £ v {z) of the universe as 
a function of frequency v and redshift z. 

The SS98 calculation of stellar emissivity closely follows this elegant analysis, with minor modifications. SS98 
also obtained metallicity correction factors for stellar radiation at various wavelengths. Decreased metallicity at high 
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redshifts gives a bluer stellar population spectrum Q , Q . 

The stellar emissivity in the universe is found to peak at 1 < z < 2, dropping off steeply at lower redshifts and more 
slowly at higher redshifts. Indeed, observational data from the Hubble Deep Field to show that metal production has 
a similar redshift distribution, such production being a direct measure of the star formation rate (see, e.g., Ref. pl|). 

With the co- moving energy density u u (z) evaluated ^] (SS98), the optical depth for 7-rays owing to electron- 
positron pair production interactions with photons of the stellar radiation background can be determined from the 
expression |lj 



t(E , z e ) 



dz ■ 



dt 
dz 



dx ■ 



dv (1 + zf 



u v {z) 
hv 



(s) 



(2) 



where s = lE^hvxil + z), Eq is the observed 7-ray energy at redshift zero, v is the frequency at redshift z, z e is 
the redshift of the 7-ray source, x = (1 — cos6>), and the pair production cross section o~ n is zero for center-of-mass 
energy yfs < 2m e c 2 , m e being the electron mass. Above this threshold, 
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where (3 = (1 — 4m 2 ! c 4 /s) 
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FIG. 5. The opacity r of the universal soft photon background to 7-rays as a function of 7-ray energy and source redshift 
(from SS98) §). These curves are calculated with and without a metallicity correction. 

Figure || shows the opacity t(Eq, z) for the energy range 10 to 500 GeV, calculated by SD98 both with and without 
a metallicity correction. Extinction of 7-rays is negligible below 10 GeV. The weak redshift dependence of the opacity 
at the higher redshifts as shown in Figure [5] indicates that the opacity is not very sensitive to the initial epoch 
of galaxy formation, z max . In fact, the uncertainty in the metallicity correction (see Figure ^) would obscure any 
dependence on z max even further. 



IV. THE EFFECT OF ABSORPTION ON THE SPECTRA OF BLAZARS AND GAMMA-RAY BURSTS 

With the 7-ray opacity t(Eq, z) calculated out to z = 3, the cutoffs in blazar 7-ray spectra caused by extragalactic 
pair production interactions with stellar photons can be predicted. Figure ^, based on the results given in Ref. j^] 
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(SS98), shows the expected effect of the intergalactic radiation grazar and 7-ray burst spectra. This figure plots the 
critical energy for absorption (i.e., for r = 1) versus redshift. For energies much above the critical energy, the optical 
depth is greater than 1, leading to a predicted cutoff in the spectrum of the extragalactic source. 




FIG. 6. The critcal energy for 7-ray absorption above which the optical depth is predicted to be greater than 1 as a function 
of the redshift of the source (from the results of SS98) Q (see text). 

The discovery of optical and X-ray afterglows of 7-ray bursts and the identification of host galaxies with measured 
redshifts (see, e.g., Refs. |35| and pq| ) has lead the accumulation of evidence that these bursts are highly relativistic 
fireballs originating at cosmological distances |57| and may be associated primarily with early star formation |38| . 

As indicated in Figure || 7-rays above an energy of ~ 15 GeV will be attenuated if they at emitted at a redshift of 
~ 3. On 17 February 1994, the EGRET telescope observed a 7-ray burst which contained a photon of energy ~ 20 
GeV p9[ |. As an example, if one adopts the opacity results which include the metallicity correction, the highest energy 
photon in this burst would be constrained probably to have originated at a redshift less than ~2. Future detectors 
such as GLAST Rcf. Q], may be able to place better redshift constraints on bursts observed at higher energies. Such 
constraints may further help to identify the host galaxies of 7-ray bursts. 

Observed cutoffs in grazar spectra may be intrinsic cutoffs in 7-ray production in the source, or may be caused by 
intrinsic 7-ray absorption within the source itself. In fact, models of quasar emission can predict natural cutoffs in 
quasar emission spectra in the relevant energy range above ~ 10 GeV. Whether or not cutoffs in grazar spectra are 
primarily caused by intergalactic absorption can be determined by observing whether the grazar cutoff energies have 
the type of redshift dependence predicted here. 
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